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We measure the hadronic contribution to the vacuum polarisation tensor, and use it to estimate the hadronic 
contribution to (g — 2) M , the muon anomalous magnetic moment. 



1. Introduction 

The vacuum polarisation Tl(q 2 ) is defined by 
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J ^ is the electromagnetic current, Q 2 

The 1-loop diagram is logarithmically diver- 
gent, n is additively renormalised. The value of 
n can be shifted up and down by a constant de- 
pending on scheme and scale without any physical 
effects, but the Q 2 dependence of U(Q 2 ) is phys- 
ically meaningful, and it must be independent of 
scheme or regularisation. 

The vacuum polarisation tensor enters physics 
in several important ways. It is responsible for 
the running of a em , which must be known very 
precisely for high-precision electro-magnetic cal- 
culations (for example of g — 2). Secondly, the 
optical theorem implies that the cross-section for 
e + e~ — » hadrons is proportional to the imagi- 
nary part of II at timelike q. This cross-section is 
usually given in terms of the ratio R, 



R(s) 



■^hadrons 



00 



<7 R +, 



00 ' 



(2) 



"Talk presented by P. Rakow at Lattice 2003 



There is an important set of dispersion relations 
connecting R and LI 
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2. The lattice calculation 



Our present calculation is carried out in the 
quenched approximation, using clover fermions 
and an 0(a) improved electromagnetic current. 
We have used three f3 values (6.0, 6.2 and 6.4) 
with several quark masses in each case. Full de- 
tails of the calculation can be found in 1 . 

IT can be split into two parts, a fermion line 
connected contribution Cn, and a disconnected 
contribution An, Fig.^ 
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We only calculate Cn, but the term An (which 
violates Zweig's rule) is expected to be very small. 

In Fig. [2] we compare our lattice data with con- 
tinuum perturbation theory. Firstly, we see very 
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Figure 1. Examples of fermion-line- connected 
and disconnected contributions to the vacuum po- 
larisation. 

good agreement between the two lattice sizes, 
showing that finite size effects are not a major 
problem. Perturbation theory works very well ex- 
cept at the largest Q 2 values, where discretisation 
errors start to show up, and at low Q 2 (below 
~ 2GeV 2 ). Both lattice sizes agree, so this is 
not a finite size effect. Can we understand this in 
terms of non-perturbative physics? 
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persion relations eq.©. These relations connect- 
ing R(s) to II(Q 2 ) are very stable, so even a crude 
model of R is likely to give a good estimate of II. 
The simplest model of R consists of (5-function 
contributions from the vector mesons (p, oj, 4>) 
and a flat continuum beginning at a threshold sq, 

R(s) = Y. e } ( AS ( S ~ m v) + BQ ( S - s o)) ■ (5) 
/ 

Using the dispersion relation gives 

C n (Q 2 ) = Bln[a 2 (Q 2 + «„)] - f 2 + K.(6) 

We have already measured m v and fy (the de- 
cay constant) so the weight and position of the 5- 
function are already known. B and sq are found 
by fitting. As seen in Fig. eq. © gives an excel- 
lent fit, and can be used to extrapolate our data 
to lower Q 2 . 
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Figure 3. Eq.$) compared with the lattice data. 



Figure 2. A check for finite size effects, compar- 
ing a 16 4 lattice (white) and a 32 4 lattice (black) 
at (3 = 6.0, k = 0.1345. The curves show con- 
tinuum perturbation theory. The dashed curve is 
at zero quark mass, the solid curve includes mass 
effects. 

One approach, introduced in uses the dis- 



3. The muon anomalous magnetic moment 

The anomalous magnetic moment of the muon 
can be calculated to very high order in QED (5 
loop), and measured very precisely, (g — 2)^ 
is more sensitive to high-energy physics than 
(g — 2) e , by a factor m 2 /m 2 , so it is a more 
promising place to look for signs of new physics, 
but to identify new physics we need to know the 
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conventional contributions very accurately. QED 
perturbative calculations take good account of 
muon and electron loops, but at the two-loop level 
quarks can be produced, which in turn will pro- 
duce gluons. The dominant contribution comes 
from photons with virtualities ~ rnt, which is a 
region where QCD perturbation theory will not 
work well. 

The traditional route for estimating these 
hadronic contributions is by using a dispersion 
relation involving the cross-section R(s). 
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By distorting the contour of integration this can 
be rewritten as 
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[Cn(Q 2 )-C n (0)] (9) 



The kernel in this relation is 
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For each data set we calculate the integral V . 
The integral is dominated by Q 2 ~ 3m 2 , so some 
Q 2 extrapolation is needed. 

In Fig. 01 we extrapolate the integral I* to the 
continuum limit and to the physical quark masses 
(m 2 s = m 2 , 2m 2 K — to 2 ) using the ansatz 
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I f = (Ai + A 2 a 2 ) + (B 1 + B 2 a 2 )m 2 ps 
Our values are 

I u = I d = 0.0389(21); I s = 0.0287(9) . (12) 

Adding these values, we get our final answer 
„2 4J « + I d + r 
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:w , {) - 446(23) x 10- 10 (13) 

This value agrees with the pioneering lattice cal- 
culation of [H|, a^ ad = 460(78) x 10~ 10 . However 
both lattice values lie lower than the experimental 
value 683.6(8.6) x 10 -10 g). Possibly this short- 
fall is due to the absence of two-pion states in the 
quenched calculation. 
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Figure 4. The results for the integral I? , and our 
extrapolations to the continuum limit. 

4. Conclusions 

We have seen that perturbation theory works 
above Q 2 - 2 GeV 2 . We can describe the entire 
Q 2 region very well with a dispersion relation us- 
ing a simple model of R(s) [2]- 

From the low Q 2 region of the vacuum polari- 
sation we can extract a lattice value for a 1 }, ad , the 
hadronic contribution to the muon's anomalous 
magnetic moment, which is of the right order of 
magnitude. Our estimate could be improved by 
using a larger lattice size, enabling us to reach 
lower Q 2 which would reduce uncertainties from 
extrapolation. Naturally, dynamical calculations 
would be very interesting. 
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